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Abstract 

We use nucleon-nucleon phase shifts obtained from experimental data, together with the chiral 
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the energy dependence produced by short-distance dynamics is well approximated by a two-term 
polynomial for Ti a b < 200 MeV. We also find that a quantitative description of NN dynamics is 
possible, at least in this channel, if one treats the long-distance parts of the chiral NN potential 
in perturbation theory. However, in order to achieve this we have to choose a separation scale R 
that is larger than 1.0 fm. 
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I. INTRODUCTION 



The chiral symmetry of the strong interaction places significant constraints on hadron- 
hadron interactions at low energies. Chiral perturbation theory (yPT) implements these 
constraints in a systematic fashion. (For recent reviews see Refs. |2J.) In contrast to the 
situation in the single-nucleon and mesonic sectors, low-energy nucleon-nucleon interactions 
do not vanish in the chiral limit, with experimental manifestations of the strength of the 
NN interaction at low energies being provided by the large np scattering length in the 1 So 
channel and the presence of a bound state (the deuteron) in the 3 Si channel. 

This complicates the application of perturbation-theory methods to multi-nucleon in- 
teractions. Weinberg's pioneering efforts 0] in the early nineties proposed surmounting 
this difficulty by making a ;^PT expansion for the NN potential, V, that goes into the 
Schrodinger equation, i.e. solving: 

(E-H )\i;)=V\i;), (1) 

with 

V = V® + V™ + + .... (2) 
Here the superscripts indicate the power of the (presumably) small quantities A "^" B , (p 
is the momentum of the NN collision and A x $b the scale of chiral-symmetry breaking) that is 
present in that piece of V. ;\PT calculation then reveals that = C 5^ 3 '(r) + includes 
one-pion exchange and a short-range interaction, V^ 2 ' includes higher-derivative short-range 
interactions, together with two-pion exchange diagrams constructed from the leading-order 
%PT Lagrangian, and involves the so-called "next-to-leading order" two-pion exchange, 

(2) 

where the ttttNN vertices are those from of ;\PT. In Weinberg's original paper there 
was no distinction made between the power counting for the piece of the potential that is 
operative at long distances (r ~ and that for the delta functions and their derivatives 

which represent the shorter-range mechanisms (r < l/m^) in this approach. Both are 
assumed to give a contribution to the NN potential of an order given by naive dimensional 
analysis in powers of p and m^. 

The implications of the expansion (T5]) for NN scattering data were first examined in 
detail by Ordonez et al. in their landmark 1996 paper j4|]. These authors analyzed NN 
scattering, not only in the 1 So, but also in a number of other partial waves. Improvements 
on this analysis were made in the work of Epelbaum and co-authors [5|, while peripheral 
NN partial waves were analyzed in first- and second-Born approximation in Ref. |6|. Two 
sets of authors have subsequently extended these analyses to fourth order in V [7, H, H|. 
In all these studies the expansion ([2]) yielded reasonably convergent results for NN phase 
shifts. 



Kaplan, Savage, and Wise (KSW) [10|, lll| have demonstrated that the assumption that 



naive dimensional analysis sets the size of the short- distance pieces of the potential is chirally 
inconsistent, in the sense that iterations of the leading-order V ([2]) via the Schrodinger equa- 
tion (pQ) generate divergences proportional to m\ [12J. But, at leading order in Weinberg's 
expansion, there is no counterterm to absorb this divergence. KSW proposed an alternative 
expansion, wherein the delta-function interaction Co8^ 3 '(r) was promoted to leading-order, 
but one-pion exchange retained its naive-dimensional-analysis scaling, making it sub-leading 



in this expansion. This expansion was carefully examined in the 1 So channel 10|, [ll|, [13 
and appears to converge reasonably well there. However, the KSW expansion does not con- 
verge in the 3 Si channel for momenta > 100 MeV [13j. A compromise proposal is to expand 
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one-pion exchange about its chiral limit value, in which case it is leading in the 3 Si and 



sub-leading in 1 Sq [14| . The resulting series has been shown to converge, albeit slowly. 



However, two-pion-exchange corrections are yet to be considered in light of the analysis 



of the 1 So of Refs. 12|, ll4(. The discussions in these works focused on establishing the 
correct power counting for the NN potential V in the case that the long-range potential is 
one-pion exchange, i.e. V = V^. (For other research that bears on the role of one-pion 



exchange in this channel see Refsjl5|,0,0. 



studies, including Refs. 



°>. (For 

'sjiEE 



29, 30 



H, 0, 0, 0, 0, 0, 0, 0, 0, 01 •) Recent 



have examined the impact of higher-order pieces of V 
on the predictions that the xPT approach to nuclear forces makes for phase shifts in the 1 So 
channel. 

Like Ref. [29j, this paper examines the role of two-pion-exchange corrections to V in the 
1 So channel. We seek to answer two questions. First, is there empirical evidence for these 
two-pion exchange contributions, and, concomitantly, what impact do they have on our 
understanding of the short-range dynamics in the NN system? Second, is there any sense 
in which and are "small"? I.e. one (or more) of the expansions proposed for NN 



dynamics in Refs. [1,0, ID, 14 is converging. The weakness of one-pion exchange in the 1 So 
channel — thanks to the absence of the tensor pieces of the NN potential that generate, e.g. 
deuteron binding — means that these questions are of crucial importance in the development 
of a quantitative theory of NN scattering in this channel. 

In contrast to Refs. 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
we do not compare predictions from a treatment of the NN potential (or 
amplitude) within a given chiral power counting with data. Instead we examine the extent 
to which phase-shift data obtained by the Nijmegen group [33j, |34j can be used to make 
inferences regarding the NN potentials that are at work for r < l/m n . To do this we invoke 
xPT with only nucieons and pions as explicit degrees of freedom to obtain V at distances 
r ~ llm-n [j, 0, 0]. (The xPT expansion for V converges quite well in this domain, see, e.g., 
Refs. jsl, |35|.) We then start with the Nijmegen PWA93 phase shift 33] at a fixed energy 
and use the long-range potential at a given order in x?T to integrate in to a finite distance 
R. The formalism for this treatment is presented in Sec. [Til Such an "outside in" approach, 
that begins with phase shifts, and uses a %PT potential to deduce information about short- 
distance dynamics in the 1 S channel, was already pursued in Refs. [l5, 22, 23, 24|. However, 



in these works only the leading-order %PT potential was employed for this purpose. 
And while the NNLO chiral potential V {0) + V {2) + was used to do the integration in 
Ref. [29|, only values R <C l/m n were discussed there. 

When such an analysis is performed over the range < T] a b < 200 MeV it yields infor- 
mation on the energy dependence that the short-range interaction must have if it is to be 
used in concert with the long-range V of %PT (at a fixed order) to reproduce data. (This 
method is very similar to those of the phase-shift analysis itself, although in Ref. [33J a 
different long-range potential was employed.) Since our effective field theory (EFT) does 
not contain explicit Delta degrees of freedom we limit ourselves to data at fairly low ener- 
gies: T lab < 200 MeV. The results for the inferred energy dependence due to short-range 
iViV dynamics are presented in Sec. IIHI While any energy dependence due to short-range 
dynamics is in principle possible, short-range potentials with coefficients that are natural 
with respect to the high scale 1/R will lead to smooth, not rapid, energy dependence on the 
interval < Tj a b < 200 MeV. As one might expect, if R is too large, the energy dependence 
over this entire energy interval is not smooth, since the NN collision probes details of the 

regulator if p cm = J ^ lab is greater than or of order -|. However, we find that for R < 1.8 
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fm only smooth energy dependence of the short-range potential is needed, and that this 
conclusion holds essentially irrespective of the ;^PT long-range potential that governs NN 
dynamics in the region from r = R to r = oo. 

The second question is whether dynamics in that region can be understood in pertur- 
bation theory. In Section IIVI we show that the long-distance pion-exchange interactions 
derived from %PT can indeed be treated in perturbation theory in the region r > 1.0 fm. 
This result is obtained via a perturbative treatment of the chiral potential used in the 
Schrodinger equation in Sec. [Til Our perturbation theory preserves the asymptotic wave 
function, and hence the phase shifts. In this way we can examine the extent to which the 
energy dependence displayed in Sec. UTTJ can be understood using perturbation theory. (A 
similar long-distance perturbation theory was developed in Ref. [29], but there it was used 
to integrate the Schrodinger equation to R ~ 0, whereas here our focus will be on finite R.) 

In Sec. |V]we present one concrete realization of short- distance physics. We adopt a two- 
parameter form for the potential in the region r < R, and obtain values for the relevant 
parameters that generate the various different energy-dependencies displayed in Sec. IIHI 
(We stress, however, that the results of Sec. IIHI for the boundary condition at r = R do not 
depend on a particular model of the physics that is operative for r < R.) Finally, in Sec. |VT] 
we offer some conclusions and an outlook. 



II. SOLVING THE SCHRODINGER EQUATION 

We choose to work in co-ordinate space and hence the basic task is to solve the radial 
Schrodinger equation: 

1 d f 2 dR{r)\ + M 



r>*V dr j+ W (E-V(r))R(r) = 0. (3) 



Here, E is the cm. energy of the NN system, and 4f its reduced mass. We adopt M 
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938.918 MeV. By the substitution u(r) = R(r)/r the equation reduces to 

d 2 u k {r) 



dr 2 



+ (k 2 ~U(r))u k (r)=0, (4) 



where k 2 = and U(r) = MV J;' - . Our objective in this section is to solve this equation 
in the region R < r. Because of this, an important feature of our solution is that we do 
not use the usual boundary condition u(0) = 0. Instead, we invert the problem and use the 
Nijmegen PWA93 [34| extraction of the 1 So phase shifts as input that provides the boundary 
condition at r — > oo: 

ship + 6(h)) 

u k [r r ^oo = — — r , 5 

5(k) being the experimental phase-shift corresponding to c.m.-frame relative momentum k. 
The normalization is done in such a manner that the asymptotic zero-energy solution is 
given by 1 — r/a, a being the scattering length. The purpose of solving the Schrodinger 
equation in this way is to obtain an energy- dependent matching condition at any R (which 
is naively of the order of or smaller than the range of pion-exchange interaction). This 
matching condition is defined by the logarithmic derivative: 

l(k-R)=\^f\] . (6) 

lU k (r)\r=R 
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For example, in the absence of any long-range potential (i.e. the case V = 0) we have: 



-f(k;R) = kcot(kR + S(k)). (7) 

If we incorporate any long-distance potential V in the analysis for r > R then the form ([7]) 
will no longer apply, and in general, we will have only numerical results for j(k; R). 

Regardless, the energy dependence of 7(/c; R) can be fitted with a polynomial in k 2 as- 

1 (k;R) = J2MR)(k 2 y. (8) 

i 

The fact that we represent 7(/c; R) in the form (JBJ) is guided by the form of the Lagrangian for 
an EFT of the NN interaction. Because of parity and time-reversal symmetries, only even 
powers of k can occur in NN contact interactions. Such local interactions (in the sense of 
quantum field theory), when smeared over a length scale R, will result in energy dependence 
that is smooth with respect to the scale R, and so we expect that Ai(R) ~ R 2,1 ^ 1 . A lucid 
description of this philosophy of regulating the short-range part of a two-body potential can 



be found in 36 



The coefficients in Eq. ([S]) are then manifestations of the short-range NN interaction in 
the channel. If one wishes these parameters can in turn be used to construct a model 
that is a particular realization of that short- distance physics. In other words, once a 
potential of a given order had been used to obtain information on 7(/c;i£), we would use 
that information to extract the coefficients of NN contact interactions <2q, C 2 , C4, etc. This 
is very much in the spirit of the Nijmegen PWA93, where Stoks et al. [33j employed square 
wells for r < 1.8 fm to produce the necessary energy dependence due to short-distance 
physics. 

But, regardless of what potential generates the coefficients Ai in Eq. ([H]), the long-range 
potential (V(r) for r > R) should produce the faster energy dependence in the phase shifts. 
Removing this more rapid energy dependence by integrating the wave function from r = 00 
in to r = R and examining j(k; R) allows us to look at the behavior of observables with 
energy that is generated by short- distance dynamics. So, although a polynomial form like 
([8]) is preferred from the EFT point of view, we are, at first, agnostic about the form of 
7(/c; R), and merely report the results obtained from the integrating-in exercise when 
potentials of different orders are adopted in the region r > R. But, before we do that, we 
first give the details of the potentials that govern dynamics in that region. 



A. The Potential 

In our analysis we adopt a xPT expansion for V(r). Since we will solve the Schrodinger 
equation in co-ordinate space, we need a co-ordinate-space representation of the correspond- 
ing pion-exchange potentials. We have adapted the expressions from Ref. [f|. The potential 
corresponding to one-pion exchange is the leading-order part of V, V^°\ and can be expressed 
as _ _ 

y<°>(r) = (ri ■ f 2 )[wf°(r)0?i ■ a 2 ) + (r)S 12 (f)], (9) 

where Si 2 (r) = 3ai • m 2 • f — a\ ■ a 2 , and 

W's"\r) = f|§^ , (!0) 
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J v ' 48vrM 2 r 3 v " ' ' ^ ' 

with x = m T r and gf^jy the ttNN coupling constant. Numerically we choose m n = 
134.98 MeV and g vN = 13.1 Q. 

Similarly, the two-pion exchange potentials which yield V^ 2 -* and W 3 ) are 

V^(r) = Re(y s (r)(ffi • B 2 ) + (n ■ r 2 )W c (r) + V T (r)S 12 (r)), (12) 

for the "leading-order" pieces of two-pion exchange, which are constructed solely out of 
vertices in C^, and 

V®(r) = Re(V c (r) + (f a ■ f 2 )W ? | 27r) (r)(5 ! 1 ■ a 2 ) + (n ■ T 2 )W^\r)S 12 (r)) (13) 

for "sub-leading" two-pion exchange, which involves contributions from C n ^. 

The different coefficient functions Vx{r) and Wx(r) — with X — C, S, T referring to the 
central, spin-spin and tensor components of the potential — are 0: 

Mr) = 12 9 8 l? f l rA {-!2xK (2x) - (15 + Ax 2 ) K x {2x)} , (14) 
W c (r) = 128 ^} 4r 4 { I 1 + 2 ^(5 + 2x 2 ) - (23 + 12x 2 )] K 1 (2x) 

+x[l + 10g 2 A -g A (23 + 4x 2 )] K (2x)\, (15) 

Vs(r) = {3x K (2x) + (3 + 2x 2 ) K^x)} , (16) 

Wt\r) = ^Tjf^ W 1 + ^)( 3 + 3 * + 2 * 2 )} > ( 17 ) 

o 2 -2a; 

^( r ) = on 2^4 6 {2cix 2 (l + x) 2 + +c 3 (6 + 12a; + 10x 2 + 4x 3 + x 4 ) } , (18) 
W?°(r) = ^^^{^(^^(S + Sx + x 2 )}. (19) 

Here, and in the calculations we report on subsequently, we have omitted the 1/M pieces 
of the potentials derived in Ref. [6J. We find that these contributions have a negligible effect 
on the results for R > 1.4 fm. This supports the power counting adopted in Refs. jl, 0], 
where M ~ A 2 , leading to a sub-dominant role for 1/M corrections in the iViV potential. 
As R is lowered below 1.0 fm, including 1/M corrections to V leads to marked differences 
with our results, which raises questions about the M ~ A 2 power counting in that domain. 
However, for R = 1.0 fm the largest change induced by the addition of 1/M corrections to 
V is only 15%. Thus, even there, they can be regarded as higher order than V^ 3 -*. 

In what follows we have used g A = 1-29 to be consistent with our choice of g n N in the LO 
potential, and f v = 92.4 MeV. For c\, c 3 and c 4 we have used two sets of values corresponding 
to the low-energy extraction of Rentmeester et. al. 31] and Entem et. al. The two sets 
are tabulated in Table [H 

We now need to project out the potentials that act in the 1 So channel. This removes the 
tensor components, leaving us with 

V^\r) = -3W^\r), (20) 
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TABLE I: Different values of the low-energy constants Cj, i = 1,3,4 (in GeV 1 ). 



ci c 3 c 4 



Rentmeester et. al. 
Entem et. al. 



-0.76 -5.08 4.70 
-0.81 -3.4 3.40 



which is our leading-order (LO) chiral potential. Similarly, 

V (2 \r) = Re(W c (r) - 3V s (r)), 



(21) 



and V^°'(r) + V^ 2 '(r) is our next-to-leading order (NLO) chiral potential. And finally, 



yielding V<® (r) + V & (r) + V^ 3 '(r) as our next-to-next-to-leading order (NNLO) chiral po- 
tential. 

III. THE ENERGY-DEPENDENT MATCHING CONDITION 

To examine how the pieces of the long-range potential manifest themselves, we 

solve the Schrodinger equation (gD using the LO, NLO, and NNLO potentials f[2"0l-f[22l. 
The output of this procedure is the energy- dependent logarithmic derivative "f(k;R), at 
different distances R. 

The results for j(k; R) are displayed as a function of lab. energy Tj a b = 2 K]f for four 
different values of R in Fig. [U As mentioned before, we present results for < Ti a b < 
200 MeV. The four different panels correspond to choosing the matching point R to be 1.0 
fm, 1.4 fm, 1.8 fm, and 2.2 fm respectively. Each panel itself contains four different curves, 
corresponding to different choices for the potential V in the region r > R. 

Let us first analyze the energy dependence of 7(fc;i?). Fig. [1] shows that '-/(k; R) has 
a smooth behavior as Tj a b increases from to 100 MeV. For all i?'s examined it could 
be represented as a polynomial in k 2 in this range. However, at slightly higher energies, 
^lab = 140 MeV, there is a singularity in the energy dependence for R = 2.2 fm. This is, in 
fact, an artifact of the singularity of the cotangent derived in Eq. (C3). The singularity occurs 
because kR « ir at these energies, which means that the lab. energy is high enough that 
we are probing details of the dynamics at r < R. It is therefore not a surprise that ^(k; R) 
cannot be represented by a polynomial for Ti a b > 140 MeV. A similar comment applies to 
the results for R = 1.8 fm, although there the singularity is at 71 a b a little above 200 MeV. 

A smooth, i.e. flatter, energy dependence is highly desirable because our intention is 
to use the energy dependence of j(k; R) to extract information about the short-distance 
physics per Eq. ([H]). If we choose R = 2.2 fm or even R = 1.8 fm, we become very limited 
in the energy-range that can be used for extraction of the short- distance physics. 

Furthermore, the lower panels in Fig. [T] make it clear that the long-range potentials are 
not having a significant impact on the energy dependence of ^(k; R) for R > 1.8 fm. In other 
words, no matter what order chiral potential we use, or even if we use no chiral potential at 
all, the result for 7(fc; R) at R > 1.8 fm looks essentially the same, although there is some 
effect due to V^°\ i.e. one-pion exchange, in j(k; 1.8 fm) at higher energies. 



V^ 3 \r) = Re(V c (r) - 3W^ n \r)) 



(22) 



7 




FIG. 1: j(k;R) (MeV) plotted against T lab (MeV) for different values of R. The dotted (black) 
curve corresponds to V(r) = 0, dashed (red) represents V(r) = V^°\r), solid (green) represents 
V(r) = y(°)(r) + vW(r) and the dot-dashed (blue) represents V(r) = V<®(r) + V^ 2 \r) + V^\r). 



Both of these phenomena can be seen in one plot if we form the dimensionless quantity 
7-R, and consider the result as a function of the dimensionless quantity kR, as well as 
dimensionless numbers formed out of scales present in the NN potential: 

lR = jR(kR,mi r R,A NN R,...), (23) 



where Ann = ^rrf HJ II, 37 is a scale that sets the strength of the iViV potential at leading 

9 a 1 " 



order, and the dots indicate the other scales that will appear if or is employed in 
the extraction of 7-R from data. In Fig. [2] we have plotted ^R against the rescaled wave 
number kR. The dotted curves in the figure are generated with R = 3.0 fm, and we see that 
they are all close to the black-dotted curve, which encodes (JT]) at this radius, and so is what 
is obtained if V = 0. In other words, at R = 3.0 fm essentially all of the 1 So phase shift is 
generated by short-range (r < 3.0 fm) dynamics. This is hardly a surprise given that OPE 
is the longest-range part of the iViV force, and m^R = 2.1 at this R. But Fig. [2] also shows 
that as we decrease the scale R — the scale that defines the demarcation between "long-" 
and "short-" distance — there is increased separation between the curve of Eq. (J7J) and the 
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0.25 0.5 0.75 1 1.25 1.5 1.75 2 

kR 

FIG. 2: Here we have plotted jR vs kR, both dimensionless quantities, for different Rs. The dotted 
curves are for R = 3.0 fm, the dashed curves for R = 1.8 fm and the solid curves for R = 1.0 fm. 
In each case the black curve is the result with V = 0, the red curve corresponds to V = the 
green curve to V = V(°) + V {2) , and the blue curve is what is obtained with V = V (°) + V ^ + . 



curves obtained when V^°\ + V^ 2 \ and + + are used for the analysis. 
At a separation scale of R = 1.8 fm pionic effects generate a larger portion of the overall 
phase-shift — although the analysis also shows that this effect comes mostly from one-pion 
exchange. But at R = 1.0 fm we see large effects from pionic dynamics. At this scale we also 
see significant differences in the results for ~/R, depending on what long-range potential is 
used for the "renormalization-group evolution" from i? = ooto-R=1.0 fm. (The formalism 
for such a renormalizati on-g roup analysis of the phase shift with respect to the scale R was 



laid out in Refs. 122, 23, 2 



Indeed, at the separation scale R = 1.0 fm, the impact of pion dynamics in the iViV 
potential on ^R is largest for the NNLO V — which includes so-called "sub-leading" two- 
pion exchange — and the result obtained for j(k; R) does not appear to be perturbatively 
close to the (black) no-long-range-potential curve. The departure from this no-long-range- 
potential result is still visible if one adopts either or + as the long-range 
potential, but it is nowhere near as large. This casts doubt on whether the hierarchy of 
effects at this separation scale is still that predicted by x?T. The issue here, of course, is 
that both and V^ 3 ) have stronger singularities as r — > than does V^°\ Thus, as r gets 
smaller they eventually dominate the evolution of 7(^5 R) with R and k. The power 
counting is only a reliable guide to the relative size of contributions provided R ~ l/m n . 

For these reasons we believe that the separation of long- and short- distance physics in 
the 1 So channel within a framework will be most effective if 1.0 fm < R < 1.8 fm. If 
R > 1.8 fm there is very little impact from chiral dynamics on j[k]R). Correspondingly, 
conversion from the variable kR back to T lab produces rapid curvature at disturbingly low 
values of the laboratory energy. At R > 1.8 fm the computation of 7 is sensitive to details 



9 



of the short-distance dynamics (and not just a few coefficients in the expansion (jSJ)) for Ti a b 
only a little above 100 MeV. In contrast, at R < 1.0 fm the short- distance dynamics has less 
role in generating the energy dependence of the 1 So phase shift; but the massive attraction 
generated by V^ 3 -* causes problems of its own, as witnessed by the very large value of jR at 
k = in the presence of "sub-leading two-pion exchange" that is inferred at this separation 
scale. 

We now take the results of Figs. [T] and [2] and extract the coefficients A and A\ in the 
expansion (jSj) of 7 (A;; R) for different choices of the long-range potential V. These results 
will be used in Section [V] when we attempt to implement specific short- distance potentials 
that reproduce the low-energy behavior of j(k; R). The results for 7(^5 R) obtained at two 
different separation scales R in the range 1.0 fm < R < 1.8 fm are presented in Table HT1 
(For the results tabulated we have have used four terms in the expansion (jSJ) to convince 
ourselves that the coefficients A 2 , A3, etc. are indeed much smaller compared to either A 
otAi. For example, A 2 ~ 10~ 3 fm 3 and A s ~ 10 -5 fm 5 for R = 1.4 fm.) The numbers in 
Table HT1 reinforce the conclusion that the next-to-leading-order two-pion exchange potential 
has the largest effect, with a pronounced impact on A Q and A\ for R — 1.0 fm. We note 
that A 2 and A3 are small, and the values of Aq and A\ shown are natural with respect to 
the scale R and the underlying scale of x?T, A x sb- 

TABLE II: The coefficients Aq (in fm -1 ) and A\ (in fm) for R = 1.0 fm and R = 1.4 fm for different 
choices of long-range potential. 



Potential 


R=1.0 fm 
A Ai 


R=1.4 fm 

A A x 


No LR pot 

y(o) 

F(0) + \/(2) 
y(0) + y{2) + y(3) 


0.042 0.328 
0.192 0.008 
0.403 -0.044 
3.101 -0.648 


0.045 -0.284 
0.130 -0.37 
0.178 -0.39 
0.496 -0.518 



Finally, if sub-leading two-pion-exchange effects are so critical one must be cautious about 
the set of low-energy constants chosen for their evaluation. Fig. [3] shows the effect of using 
different values of these q's in our calculation. The two choices are those tabulated in Table [H 
It is obvious that the Nijmegen set of q's produces a stronger effect in 7(^5 R), which is not 
surprising given that the combination C3 — 2C4 sets the size of the leading singularity in 
in this channel, and this combination is much larger for the q's adopted in Ref. [3lj. It 
is interesting to note that both choices lead to a similar shape for j(k;R), although the 
magnitude is larger in the case of the Rentmeester et al. choice. Presumably the fact that 
7(fc; R) is smooth in both panels allows either choice of q's to yield a reasonable fit for 
NN scattering data that is sensitive to the 1 So phase shift. It is just that the separation 
between the long-distance dynamics in and the short-distance dynamics encoded in jR 
will differ, depending on the value of the q's that is adopted. From this point on all of the 
results that we present use the Nijmegen set of q's. 

IV. ARE PION EXCHANGES PERTURB ATIVE? 

In the previous section we showed that the logarithmic derivative 7 (A:; R) is a logical way 
to look at the short- distance parts of the iVW interaction. In this section we turn to the 
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FIG. 3: Effect of different sets of q's from Table [J on the energy-dependent matching condition at 
R= 1.0 fm and 1.4 fm with V(r) = V^°\r) + V^{r) + V^\r). 

question of whether the long-range part of V can be treated perturbatively. To do this we 
analyze how close a perturbative treatment of the potential is to the full solution of 
the Schrodinger equation in the region r > R. 

Our objective is still to solve the Schrodinger equation: 



d 2 u(r) 
dr 2 



+ (k 2 - U{r))u{r) = 0, 



(24) 



in the region r > R. Here U(r) = is obtained from the potential in this region of 

space, which in our case is the LO, NLO or the NNLO potentials generated in ;\TT. In order 
to build up the solution, let us first rewrite Eq. ff24l as: 



d 2 u k (r) 
dr 2 



+ aU(r)u k (r) = k 2 u k (r), 



(25) 



where a is a parameter that is used to describe the perturbative expansion in powers of 
V (note that this is not the same as the expansion for V itself that was described in 
Sec. H]), and the subscript k refers to the wave number of interest. The full solution can now 
be written as a power series in a as: 



Uk {r) = Y,a n u ( k \r) = uf{r)+au k l \r) + a 2 u k 2 \r) + 



,(o) 



C 1 )/ 



.2„.( 2 ), 



(26) 



n=0 



where, vS°\r) is the zeroth-order solution (in the absence of V(r)), u^'(r) is the first-order- 
in- V correction, u^ 2 \r) is the second-order correction, and so on. 

Substituting Eq. (|26l) into Eq. (|25|) and equating powers of a we get: 



dr 2 



u {0 \r) + k 2 u { k \r) = 0, 



dr- 



■ u[ n+1 \r) + k 2 u { k n+1 \r) = U(r)u[ n \r)- n > 0. 



(27) 
(28) 
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The solution of Eq. (T27j) that reproduces the PWA93 1 Sq phase shift S(k) is: 



sin (Ax + 5(h)) 
sin S(k) 



(29) 



Equations (12911 and ( 1281) can now be used, together with the Green's function technique, to 
calculate the corrections to u at first order and second order in perturbation theory. The 
Green's function for Eq. (1251) that preserves the form (1291) as r —>■ oo in the full solution 
u k (r) is: 

sin k(r-r') -r „ / 

it r < r , 

if r > r'. 



G{r,r';k) 



k 




(30) 



Therefore the first-order correction to the wavefunction can be expressed as: 



M 
J? 



G(r, r';k)U(r')uf\r')dr' 

sin Hr — r') . sm(kr' + S(k)) 
—V[r)- 



dr'. 



k ' sin 5(k) 

The second-order correction to the wavefunction is then calculable as 



uf\r) = (— ) / G(r,r';A0VV)4V)dr'. 



(31) 



(32) 



Comparison of results from Eqs. ( 1321) and ( 13T1) will allow us to assess the convergence of our 




"long-distance perturbation theory" . This perturbation theory was set up in Refs. |16|, 
and was used to discuss the convergence of the chiral expansion. In contrast to Ref. 
which integrated Eqs. (13T]) and (1321) to R = 0, we will only integrate them to some finite R. 
Hence, we once again examine the logarithmic derivative at a radius R, only this time we 
define a version that can be computed using long-distance perturbation theory: 



^ n \k-R) 



iu 



(o) 



+ 



u 



(») 



)'{r) 



u 



(0) 



+ 4" 



+ 



u 



(n) 



(33) 



r=_R 



Strictly speaking this 7^ includes effects of V which are of an order higher than n, but it 
has the advantage that it is straightforward to compute. It also results in smooth energy 
dependence when the limit R — 
To test the usefulness of Eq. 



24. 30 



is taken 

(133]) we choose the strongest potential, which, for distances 
R < 1.8 fm, we know to be the NNLO potential V(r) = V (0 \r) + V (2 \r) + V (3 \r), and com- 
pare the results with those of the previous section. We again adopt the q's of Rentmeester 



et al. [3l|]. The results are plotted in Fig. HI The dotted (black) curve is the zeroth-order 
solution, which corresponds to using the free Schrodinger equation for the integration from 
r = 00 to r = R. This, then, is the same as the "No-long-range potential" case of the 
previous section. Meanwhile, the dashed (red) curve is the first-order perturbation theory 
result, the solid (green) curve is the second-order result, and the dot-dashed (blue) curve is 
from the full solution of the Schrodinger equation obtained in Sec. IIHI 

The right panel of Fig. H] shows that — as one might expect given the results of that 
section — at R = 2.2 fm this potential has almost no impact on 7 (A;; R), and the result for 
any finite order of perturbation theory is the same. The left panel of Fig. H] reiterates that 
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FIG. 4: 7(/c; i?)(MeV) vs. T; a b(MeV) calculated in perturbation theory and compared to the full 
solution. The dotted (black) curve is the zeroth order solution, the dashed (red) curve is the 
first order perturbation theory result, the solid (green) curve is the second order result and the 
dot-dashed (blue) curve is from the full solution of the Schrodinger equation. 

at R = 1.4 fm the 7(fc; R) results with the NNLO potential are very different from the "No- 
long-range potential case" . But it also shows that this difference is almost entirely due to the 
first-order perturbation-theory correction. There is a 10% difference between the first-order 
result, 7^ and the full 7 at ~ 0, but this 10% difference has almost completely disappeared 
by 21ab = 100 MeV: perturbation theory works better at higher energies. And even at the 
lower energies where there is some noticeable deviation, the inclusion of the second-order 
correction brings us to within better than 1% of the full result for 7(^5 R). This success of 
perturbation theory implies that the coefficients Aq and A\ defined in Eq. ([8]) can instead 
be extracted using the perturbative formulae (1291 . (I3TI) . and (I32I) — rather than from the full 
Schrodinger equation solution — without introducing undue error. 

TABLE III: The coefficients Aq (in fm -1 ) and A\ (in fm) for R = 1.4 fm extracted using Eq. (|33p . 
The second and third (fourth and fifth) columns show the Aq and A\ corresponding to 7W (7^)- 
The sixth and seventh columns repeat the result of Table IIIIl for comparison. 



Potential 


7 (i) 
A A 1 


7 (2) 

A A 1 


Full soln. 

A Q Ax 


y-(o) + y-(2) 

V (0) + y(2) + y(3) 


0.130 -0.369 
0.18 -0.401 
0.569 -0.624 


0.130 -0.369 
0.178 -0.391 
0.516 -0.534 


0.130 -0.37 
0.178 -0.39 
0.496 -0.518 



In Table IIIIl we do exactly that, using the perturbative formula (133]) to determine A 
and Ai for the case R — 1.4 fm. The numbers agree well with those in Table [Til for the 
"full solution". A first-order-perturbation-theory extraction is sufficient for V(r) = 
and V(r) = + V^ 2 \ For V(r) = + + a second-order-perturbation-theory 
calculation does a noticeably better job. 
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Since the potentials get stronger as r — > the convergence of perturbation theory is 
rather slow for R = 1.0 fm. We find that for R = 1.0 fm, 7^^(fc;i?) overshoots , y(k;R) by 
almost a factor of two at k = 0. The second-order correction brings ^ 2 \k; R) to within 10% 
of the "exact" result at k — 0, and the agreement is better than 5% at T lab = 100 MeV. 
So perturbation theory formally converges for R = 1.0 fm, but the "chiral" nature of the 
perturbation theory is very much in question, since — as emphasized in Sec. IHH — the effects 
of are much larger than those of in this region. 

We conclude that the impact of the pion-exchange potentials generated in can be 
calculated in perturbation theory, provided that the region of r where that perturbation 
theory is applied is chosen judiciously. It is possible that a more sophisticated perturbation 
theory works even if R < 1.0 fm, but what is already clear from these results is that standard 
"long-distance" perturbation theory is applicable and useful in the domain R > 1.0 fm. 



V. SHORT-DISTANCE PARAMETERS AND POTENTIAL WELLS 

We now change gears and discuss how the information gleaned from the energy-dependent 
matching condition 7(/c; R) can be used to extract useful results regarding the short- distance 
physics that is operative at r < R. 

The coefficients defined in Eq. (jSj) and listed in Table [Til for different choices of R and 
V(r) are numbers. But, they can be used to build a short-distance potential that is, in effect, 
a short- distance regulator for our NN interaction. Our results suggest that two terms are 
sufficient to ensure convergence in the expansion in Eq. (jSj). Thus, any short-range potential 
we build having two parameters can be correlated to the coefficients A and A\. We have 
chosen to design our short-range potential to be a well with a bottom that has a slope so 
that 

V SR (r) = -V + fir, for r < R. (34) 

Here, Vb (the strength of the short- distance potential at r = 0) and the slope \i are the two 
parameters of our regulator which we denote as short- distance parameters (SDPs). In order 
to extract these two parameters we solve the Schrodinger equation for r < R such that the 
following condition is satisfied: 



u in {r) 



T=R 



l(k;R). 



(35) 



The motivation in choosing the short-range potential in this form is that through an appro- 
priate change of variable from r to x, where 



x 



( M\i 



E + Vp 



(36) 



the Schrodinger equation for r < R can be reduced to the Airy equation: 

cPuix) 



dx 2 



xu(x) = 0. 



The solution to Eq. ( 1371) is a linear combination of Airy functions: 

u(x) = a\Ai(x) + a 2 Bi(x). 



(37) 



(38) 
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The Airy function Bi(x) diverges at x — > oo, but we are interested in the region < r < R, 
and so must keep both solutions. a\ and 02 are arbitrary constants that are evaluated from 
boundary conditions u(0) = and u'(0) = 1. Thus changing the variable back to r we 



obtain- 



( h 2 \ s 

w(r)=7r(^-J [-Bi(x )Ai(r + x ) + Ai(x )Bi + x )] (39) 

where, x = f = (f^r. 

From this solution at r < i? we can easily obtain its logarithmic derivative at r = R: 

u'(R) ( H 2 Y Bi(x )Ai'(y ) - Ai(x )Bi'(y ) 



u(R) \MnJ Bi(x )Ai(y ) - Ai(x )Bi(y ) 



(40) 



with y = (^) 3 [R— E+v ° ] ■ This logarithmic derivative is a function of our SDPs, Vq and /i. 
It should be equal to the logarithmic derivative ^(k; R) computed in Sec. IHII Hence, using 
the lowest two terms in the Taylor expansion of that ^y(k; R), A and Ax, we can extract Vq 
and /1. 



TABLE IV: The short-distance parameters /i (in MeV fin -1 ) and Vq (in MeV) for R = 1.4 fm with 
different choices of long-distance potential. (Note that in our convention Vq < corresponds to a 
repulsive potential.) The second and third columns give the results that match the Taylor-series 
coefficients Aq and A\ of 7(fc;i2) defined in Eq. ([6]). The fourth and fifth (sixth and seventh) 
columns provide SDPs that match the Aq and A\ for 7W (7^)) of Eq. (|33p . 



Potential 


Pull soln. 
V Vq 


7 (i) 
V- Vq 


7 (2) 

/" Vq 


No LR pot 

y(o)( r ) 

V<°)(r) + V (2 \r) 
V (°)( r ) + V^{r) + V^{r) 


9.85 59.35 
-54.5 -10.68 
-56.0 -17.39 
-65.0 -54.4 


9.85 59.35 
-54.5 -10.56 
-55.5 -18.23 
-64.0 -55.65 


9.85 59.35 
-54.5 -10.56 
-56.0 -16.94 
-65.0 -51.1 



To demonstrate, we have calculated the short- distance parameters, Vq and /z, in this way 
for R — 1.4 fm. The results are given in Table HVl Results for a different two-parameter 



short- distance potential and the case of V = and V = V^ can be found in Ref. 15 

We had commented earlier in Sec. II VI that A Q and Ai can be extracted from the matching 
condition calculated in perturbation theory. Consequently, one can then proceed to extract 
\i and Vq from these "perturbative" Aq and A\. The results are shown in the third to sixth 
column of Table [IV] and are gratifyingly close to those obtained from the full Schrodinger 
equation solution of Sec. IHII 

The results of Table IIVI show that a nice quantitative description of the iVTV interaction 
in the 1 So channel up to Tj a b = 200 MeV can be obtained by using perturbation theory for 
r > R to calculate ^(k; R), then using that information to determine the parameters present 
in a simple short-distance potential such as 
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VI. CONCLUSION 



Our analysis of NN scattering in the 1 So channel is similar in its philosophy to that 
advocated by Lepage [36], in that we separate the potential into a long-distance part- 



determined by xPT — and a short-distance part, which can be parameterized in a variety 
of ways without affecting the predictions for NN scattering. Indeed, we have gone further 
than Lepage, and argued that the energy-dependent logarithmic derivative at R, j(k;R), 
is a convenient way to summarize information about the impact of physics that is short- 
range with respect to the scale R. This information on the short-range NN dynamics can 
be obtained from the NN phase shifts if we know the long-range potential. Combining 
%PT long-range potentials (computed to different orders in the chiral expansion) with the 
experimental phase-shifts from the Nijmegen PWA93 we obtained ^(k; R) for a range of 
scales R. 

For the extraction of model-independent information on short-distance physics it is de- 
sirable that j(k; R) be a smooth function of k. This requirement limits the energy range 
within which one can operate if R is chosen to be greater than 1.8 fm. For values of R in 
the range from about 1.0 fm to 1.8 fm 7 (A;; R) can be well described by a polynomial in 
k 2 over the entire range < Tj a b < 200 MeV. The coefficients of the terms in this poly- 
nomial represent the effect of the short- distance NN interaction. We have found that for 
1.0 fm < R < 1.8 fm the first two terms of this polynomial expansion are sufficient to ensure 
a good representation of the energy dependence of 7. So, in this range of R, details of the 
short- distance potential are not important: its effects can be summarized in A Q and A]_. 

The long-ran ge p otential V determines the evolution of j(k; R) with the separation scale 
R (see Refs. [231. l24f for a derivation of the renormalization-group equations associated with 
this evolution). The scale R in our analysis plays a similar role to the scale A in the potential 
Viowfc [HI]. In that approach A is a cutoff that separates the momentum-space states explicitly 
included in the Hilbert space from those states whose effects are encoded in V. There is then 
evolution of the NN potential V with A so that S-matrix elements (which are equivalent 
to phase shifts) stay the same as the momentum-space cutoff is changed. In our work the 
short-distance part of the potential also evolves in such a way that the 1 S phase shifts 
always agree with those from the Nijmegen PWA93. This evolution can be traced over a 
range from R = 00 to R = 1.0 fm. The fact that, at the lower end of this range, details 
of the short- distance potential do not play a key role in describing NN phase shifts below 
Tub = 200 MeV, is presumably a corollary of the existence of a "universal" V\ ow k if the scale 
A is lowered to 2 fm -1 . 

In the energy range < Tj a b < 200 MeV we find that the one-pion exchange interaction 
is the dominant effect in this evolution for R > 1.8 fm, but its impact on j(k] R) is not 
sizeable. We obtain a larger shift in 7(^5 R) compared to its value in the absence of any 
long-range potential if we consider sub-leading pieces of the xPT V. In particular, if we 
compute with V^ 3 \ and consider R < 1.4 fm, we see a significant impact of pion-exchange 
dynamics on the 7 (A;; R) inferred from the NN phase shifts. 

But, even in this long as we maintain R > 1.0 fm, we find that the evolution 

of 7(fc; R) due to long-distance effects can be understood in perturbation theory. Hence we 
can use perturbation theory in the xPT potentials in the region 1.0 fm < R < 00 to connect 
experimental data to information on the short- distance dynamics in the region R < 1.0 
fm. The short-distance part then contains non-perturbative dynamics, but these effects 
are summarized in a few coefficients in the polynomial expansion of the energy- dependent 
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matching condition. Those coefficients can in turn be matched to an explicit realization of 
short- distance physics, e.g. the one discussed in Sec. [V] 

This represents an updated version of the proposal first made by Kaplan, Savage, and 
Wise in Refs. (13, [lj: that the contact interaction parameterizing short-distance NN dy- 



namics should be promoted to order Q 1 so that it is a leading-order effect, with pion 
exchanges then being sub-leading. This proposal was examined in detail for NN scattering 



in Refs. 13|, ll4|, and was found to work well in the 1 So channel as far as the pion potential 
was concerned. We find that can also be treated in perturbation theory, and its 
effects are smaller than V<® for r > 1.0 fm. The NNLO piece of V, , can also be treated 
in perturbation theory, but only if it is regulated at a scale of order, or well below, the 
chiral-symmetry breaking scale A x sb- The large size of corrections due to raises the 
question of whether the chiral expansion for V is behaving as decreed by the ;\PT power 
counting. 

In particular, some authors have suggested that an expansion in which the Delta-nucleon 
mass difference, A, is treated as a light scale, might have better convergence properties [H, 

H, Hi H H m • 

Our results support this view — at least if one takes R to be small 
enough that chiral dynamics plays a significant role. More recently Robillota has pointed 
out that the poor convergence of xPT without an explicit Delta for the scalar-isoscalar piece 
of the NN potential is related to issues with the description of the nucleon's scalar form 



factor [45|. He observes that in that case too, contributions from Delta loops exceed those 
from nucleon loops once one considers distances r of 1 fm or less. 

We close with some possible implications of our results for future partial-wave analyses 
(PWAs). Here we have used the Nijmegen-extracted phase shifts 34j and our separation 
of NN dynamics into lon g-ra nge and short-range parts is similar to that adopted in the 
Stoks et. al. 1993 PWA [33]. However, in Ref. [33j an upgraded version of the Nijm78 
potential was used to describe the long-distance dynamics. Minimally, this paper shows 
that in the 1 So channel we can use xPT potentials and still get sensible energy dependence 
in our matching condition which encodes the short-distance NN dynamics. We have shown 
that the region 1.0 < R < 1.8 fm yields best results in this regard. This is in accord with 
the more recent PWA for which details are as yet unpublished [3lL HI. But our results 



suggest that — at least in the 1 S channel — future PWAs using these values of R could employ 
perturbation theory to compute the effects of the long-distance potentials obtained in ;\RT. 
It might also be possible to lower the separation scale R below 1.0 fm and encode the smooth 
energy dependence due to dynamics at r < R in energy- dependent square wells, but then 
the parameters obtained for these wells would be very different depending on whether the 
LO or NNLO NN potential were used — or even presumably, if different sets of q's were 
used in computing V^ 3 \ Extensions of the methods laid out here to other partial waves are 
straightforward and would facilitate such a future EFT-based partial-wave analysis. 
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